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Abstract. This talk reports work done in collaboration with Jin Hur, Choonkyu Lee 
and Hyunsoo Min concerning the computation of the precise mass dependence of the 
fermion determinant for quarks in the presence of an instanton background. The result 
interpolates smoothly between the previously known chiral and heavy quark limits of 
extreme small and large mass. The computational method makes use of the fact that 
the single instanton background has radial symmetry, so that the computation can 
be reduced to a sum over partial waves of logarithms of radial determinants, each of 
which can be computed numerically in an efficient manner using a theorem of Gelfand 
and Yaglom. The bare sum over partial waves is divergent and must be regulated and 
renormalized. We use the angular momentum cutoff regularization and renormalization 
scheme. Our results provide an extension of the Gelfand- Yaglom result to higher 
dimensional separable differential operators. I also comment on the application of this 
approach to a wide variety of fluctuation determinant computations in quantum field 
theory. 
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1. Introduction 

In this talk I present an overview of some recent progress PP in computing one loop 
quantum vacuum polarization effects. Mathematically, this requires computing the 
determinant of a fluctuation operator, which describes the quadratic fluctuations about 
a semiclassical solution. This is a difficult problem, but it is worth studying as it has 
important physical applications to computations of the effective action, the partition 
function or the free energy. It is also an interesting mathematical problem to learn 
about the spectral properties of partial differential operators. In this talk I focus 
mainly on the computation of fermion determinants in nontrivial background fields, 
which is an important challenge for both continuum and lattice quantum field theory. 
Explicit analytic results are known only for very simple backgrounds, and are essentially 
all variations on the original work of Heisenberg and Euler [3 El- For applications 
in quantum chromodynamics (QCD), an important class of background gauge fields 
are instanton fields, as these minimize the Euclidean gauge action within a given 
topological sector of the gauge field. Furthermore, instanton physics has many important 
phenomenological consequences jH 13 Ej- Thus, we are led to consider the fermion 
determinant, and the associated one- loop effective action, for quarks of mass m in an 
instanton background. Here, no exact results are known for the full mass dependence, 
although several terms have been computed analytically in the small mass |3 IZl IHl and 
large mass E] limits. Recently P, with J. Hur, C. Lee and H.Min, we presented a 
new computation which is numerical, but essentially exact, that evaluates the one-loop 
effective action in a single instanton background, for any value of the quark mass (and 
for arbitrary instanton size parameter). The result is fully consistent with the known 
small and large mass limits, and interpolates smoothly between these limits. This could 
be of interest for the extrapolation of lattice results (TU] , obtained at unphysically large 
quark masses, to lower physical masses, and for various instanton-based phenomenology. 

Our computational method is simple and efficient, and can be adapted to many 
other determinant computations in which the background is sufficiently symmetric so 
that the problem can be reduced to a product of one- dimensional radial determinants. 
While this is still a very restricted set of background field configurations, it contains 
many examples of interest, the single instanton being one of the most obvious. The 
method is based on the Gelfand-Yaglom method for computing determinants of ordinary 
differential operators IT ^ IT! ^ IT ^ ITH | But in higher-dimensional problems with 
partial differential operators, it is known jT3] that the naive generalization is divergent, 
even for simple separable problems where there is an infinite number of 1-D determinants 
to deal with. Physically, this divergence reflects the fact that in dimensions greater than 
one, one must confront renormalization. Our results can be viewed, in fact, as giving 
an extension of the Gelfand-Yaglom result to higher dimensional separable differential 
operators. 
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2. Preliminaries: Effective Action in an Instanton Background 

The first step is to recall that since the instanton background field is self-dual, we can 
deduce the fermion determinant from a computation of the determinant of the associated 
Klein-Gordon operator. This is because self-dual gauge fields have the remarkable 
property that the Dirac and Klein-Gordon operators in such a background are isospectral 
jUEl- This implies that the renormalized one- loop effective action of a Dirac spinor field 
of mass m (and isospin |), T^^^{A]m), is related to the corresponding scalar effective 
action, r^^^{A]m), for a complex scalar of mass m (and isospin |) by |31IH] 

ri^{A- m) = -2 rL(A; m) - i In (^^^ , (1) 

where n is the renormalization scale. The In term in corresponds to the existence of a 
zero eigenvalue in the spectrum of the Dirac operator for a single instanton background. 
We consider an SU(2) single instanton background, in regular gauge IUCH] : 

A (x) = A^ix)— = ^^'^"^"^^ F (x) = (x)— = (2) 

where rj^^a are the standard 't Hooft symbols [H . 

The one-loop effective action must be regularized. We use Pauli-Villars 
regularization [with heavy regulator mass A], adapted to the Schwinger proper-time 
formalism, and later we relate this to dimensional regularization, as in the work of 't 
Hooft The regularized effective action has the proper-time representation 

Ti{A-m)= - r^(e-™''^-e-^'^) ( d^x ti{x\e~-<~^"'> - e- \x) , (3) 



where = D^D^, and = — iA^{x). The renormalized effective action, in the 



minimal subtraction scheme, is defined as |31 |H] 



rfen(^; m) = lim 

A^oo 



lim 

A^oo 



ri(A;H-l^ln(^)/.^xtr(F,.F^ 
Ti{A-m)- 




(4) 



where we have subtracted the charge renormalization counter-term, and p is the 
renormalization scale. By dimensional considerations, we may introduce the modified 
scalar effective action r^^(mp), which is a function of mp only, defined by 

rL(A; m) = fL(mp) + i In(pp) , (5) 

and concentrate on studying the mp dependence of F^j^ (mp). Then there is no loss of 
generality in our setting the instanton scale p = 1 henceforth. 

It is known from previous work that in the small mass j3J [71 |H] and large mass El 
limits, rfgn(m) behaves as 

{a(l/2)-h^(lnm-h7-ln2)m^-F... , m^O 
Inm 1 17 232 7916 
— — -h — H" • * * Tfl — ^ OO 
6 75m2 735m4 2835m6 148225m8 




Figure 1. Plot of the analytic small and large mass expansions for fj^jj(m), from 
Equation ©. Note the gap in the region 0.5 < m < 1, in which the two expansions 
do not match up. 

where a(l/2) = -A _ 2C'(-1) - |ln2 ~ 0.145873, and 7 ^ 0.5772... is Euler's 
constant. This small mass expansion is based on the fact that the massless propagators 
in an instanton background are known in closed-form. On the other hand, the large 
mass expansion in © can be computed in several ways. A very direct approach is to 
use the small-s behavior of the proper-time function appearing in (jnj, as given by the 
Schwinger-DeWitt expansion. 

Equation (jH)) summarizes what is known analytically about the mass dependence of 
the renormalized one-loop effective action in an instanton background. This situation 
is represented in Figure ^ which shows a distinct gap approximately in the region 
0.5 < m < 1, where the small and large mass expansions do not match up. 

3. The Gelfand-Yaglom Technique 

A beautiful result of Gelfand and Yaglom [TT] provides a spectacularly simple way to 
compute the determinant of a one-dimensional differential operator. Their result has 
since been extended in various ways [T^ ITH| ITi| IT^ ITT)] , but here we will only need 
their result for radial differential operators. Suppose Aii and are two second 
order ordinary differential operators on the interval r G [0, 00), with Dirichlet boundary 
conditions assumed. In practice we will choose A^i to be the operator of interest, and 
Ai2 to be the corresponding free operator. Then the ratio of the determinants is given 
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by 

detail ^ (M^\ 

detM2 R^[MR)J 
where ipi{r) (for i = 1,2 labelhng the two different differential operators, Aii and A^2) 
satisfies the initial value problem 

MiiJi{r) = 0, with i)i{0) = and ^-(0) = 1 • (8) 

Note that no direct information about the spectrum (either bound or continuum states, 
or phase shifts) is required in order to compute the determinant. All that is required is 
the integration of (jHl), which is straightforward to implement numerically. 

4. Radial Formulation 

Returning now to the instanton determinant problem, we can use the fact that the 
single instanton background (j21) has radial symmetry [4j. Thus, the regularized one- 
loop effective action Q can be reduced to a sum over partial waves of logarithms of 
determinants of radial ordinary differential operators. Each such radial determinant 
can be computed using the Gelfand-Yaglom result (|7j). Unfortunately, the sum over all 
partial waves is divergent. The physical challenge is to renormalize this divergent sum. 

In the instanton background with scale p = 1, the Klein-Gordon operator —D^ 
for isospin ^ particles can be cast in the radial form |1] 



92 3 d ^ 4/(/ + l) ^ 4(j-0(j + / + l) 3 



(9) 



where / = 0, |, 1, |, ■ ■ ■ , and j = | / ± | |, and there is a degeneracy factor of 
{21 + l)(2j + 1) for each partial wave characterized by (/, j)-values. In the absence of 
the instanton background, the free operator is 



92 3d 4/(/ + l) 



^^iy'2i ly ^^'y <y'2i 



(10) 



This radial decomposition means that we can express the Pauli-Villars regularized 
effective action also as 

^sr, ^ ^ ,n/i /det[7^(,+ i,,)+m2]' 

T^iA-.m) = > deg(/) < In — — /' -— + In - — ;, ^ ' 

iJ^ I V det[7^5]^ + m2] J \^ det[7i:[j^\^ +m2 



2 ' 

_ / det[7^,,,.)+A2] \ _ / det[7^,,.,, + A2] 

det[7^Jf + A2] ) '""y det[7^J-^^ + A2] ^''^ 

Here we have combined the radial determinants for (/, j = /+|) and = (/+|) — ^), 

which have the common degeneracy factor deg(/) = {21 + 1)(2Z + 2). 

The Gelfand-Yaglom technique provides a simple and efficient numerical technique 
for computing each of the radial determinants appearing in (jllj) . But to extract the 
renormalized effective action we need to be able to consider the A — > 00 limit in 
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conjunction with the infinite sum over /. This can be achieved as follows. Split the 
I sum in flllj) into two parts as : 

L oo 

TliA;m)= ^ rl(,)(A;m)+ ^ ri,(,)(A;m) (12) 

« = 0,i,... '=^+1 

where L is a large but finite integer. The first sum involves low partial wave modes, and 
the second sum involves the high partial wave modes. We consider each sum separately, 
before recombining them to obtain our final expression (j22I)- 



5. Low partial wave modes: determinants computed using Gelfand-Yaglom 

The first sum in fll2j) is finite, so the cutoff A may be safely removed, and for each / the 
determinant can be computed using the Gelfand-Yaglom result ((Tj). We can simplify 
the numerical computation further by noting that for the free massive Klein-Gordon 
partial- wave operator, 7i[;)"^ + (with TC^^i^^ given in (fTU|)). the solution to ^ is 

V,-(r) = ^H±^ . (13) 

This solution grows exponentially fast at large r, as do the numerical solutions to (jHI) 
for the operators Ti-^ij) + m?, with 'H{i^j) specified in (0). Thus, it is numerically better 
to consider the ODE satisfied by the ratio of the two functions 

^(U)(r) = ||^ ; 7^,,)(0) = l ; 7^'(,,,(0) = . (14) 

This quantity has a finite value in the large r limit, which is just the ratio of the 
determinants as in (jZj). In fact, since we are ultimately interested in the logarithm of 
the determinant, it is more convenient (and more stable numerically) to consider the 
logarithm of the ratio 

%j)(r) = ln7^(^J)(r) , (15) 

which satisfies the differential equation 



dr^ \ dr I \r hi+i{'mr) J dr ' 



with boundary conditions 

%,)(r = 0) = , S(',,,.)(r = 0) = . (17) 
The 'potential' term U(i^j){r) in (fTB|) is given by 

(18) 

To illustrate the computational method, in Figure El we plot S^i^ij^^-^{r-) and 
>S'(^_,_i ;)(r) for various values of /, with mass value m = 1 (which is in the region in which 
neither the large nor small mass expansions is accurate). Note that the curves quickly 
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Figure 2. Plots of the r dependence of 5'(; ;^i-)(r) and S'(;^i ;-)(r), solutions of the 
nonlinear differential equation in 1)16(1 . for m = 1, and for I = 0, 10, 20, 30. The upper 
curves are for S(^i ;_|_i')(r), while the lower ones are for 5'(;+i 
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Figure 3. Plot of the I dependence of P{1) = S'^; — oo) + 5'(;+^ i){r = oo), for 
TO = 1. P{1) behaves like 0( j) for large I. 



reach an asymptotic large-r constant value, and also notice that the contributions from 
S(^l^l_^_^){r = oo) and S'^^^i ;-)(r = oo) almost cancel one another when summed. Indeed, 

^(0 = '5(M+i)(^ = oo) + 5(;+i_i)(r = oo) 

~oQ) , /-.oo (19) 
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This behavior is illustrated in Figure El At first sight, this looks like bad news, because 
when including the degeneracy factor the sum over I is 

Y: Tf{A;m)= (2/ + 1)(2/ + 2)P(0 (20) 

Z=0,i,- i=0,i,... 

Therefore, the sum ()2()|1 has potentially divergent terms going like L^, L and InL, as 
well as terms finite and vanishing for large L. Remarkably, we find below that these 
divergent terms are exactly canceled by the divergent large L terms found in the 
previous section for the second sum in (fT^ . 



6. High partial wave modes: determinants computed using radial WKB 



In the second sum in (fT^ we cannot take the large L and large A limits blindly, as each 
leads to a divergence. To treat these high / modes, we use radial WKB [19 , which is a 
good approximation in precisely this limit. This means we can compute analytically the 
large A and large L divergences of the second sum in (fT^ , using the WKB approximation 
for the corresponding determinants. This is a straightforward computation, the details 
of which can be found in ^] . We find the following analytic expression for the large 
L behavior 



E 

l = L+: 



i A,(Z) 



{A;m) 



1 



1 



m 



-lnA + 2L^ + 4L- - + — InL 



6 



+ 



6 



127 



1 , 

- In 2 + — 
3 2 



m In 2 H In m 



(21) 



72 3 2 2 

It is important to identify the physical role of the various terms in (12 Ij) . The first 
term is the expected logarithmic counterterm which is subtracted in (^, and explains 
the origin of the |ln/i term in (0). The next three terms give quadratic, linear and 
logarithmic divergences in L. We shall show in the next section that these divergences 
cancel corresponding divergences in the first sum in (fT^ . which were found in our 
numerical data. It is a highly nontrivial check on this WKB computation that these 
divergent terms have the correct coefficients to cancel these divergences. Note that the 
InL coefficient, and the finite term, are mass dependent. 



7. Putting it all together 

We now combine the numerical results for the low partial wave modes with the radial 
WKB results for the high partial wave modes to obtain the minimally subtracted 
renormalized effective action rfg^(m) as 





Figure 4. Plot of our numerical results for ff^^-^{m,) from (|22|l . compared with the 
analytic extreme small and large mass limits [dashed curves] from 10 . The dots denote 
numerical data points from l|22|l . and the solid line is a fit through these points. The 
agreement with the analytic small and large mass limits is very precise. 



The remarkable observation is that the large L divergences found in the numerical results 
are precisely cancelled by the analytic large L divergence found using WKB, leaving a 
finite answer for the renormalized effective action. This holds for any mass m. In 
Figure m we plot these results for r^„(m), and compare them with the analytic small 
and large mass expansions in ©. The agreement is spectacular. Thus, our expression 
(|22p provides a simple and numerically precise interpolation between the large mass and 
small mass regimes. 

As an interesting analytic check, our formula (j22p provides a very simple 
computation of 't Hooft's leading small mass result. When m = 0, the P{1) can be 
computed analytically: P{1) = In . Then 



17 1 1 

1 



+ 2L' + AL 



1, r 

- In L H 

6 72 



In 2 



a 



0.145873... 



(23) 



which agrees precisely with the leading term in the small mass limit in (jHl). 
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8. Comparison With The Derivative Expansion 

In [19], the renormahzed effective action rfg^(m) was computed using the derivative 
expansion. Recall that the philosophy of the derivative expansion is to compute the 
one-loop effective action for a covariantly constant background field, which can be done 
exactly, and then perturb around this constant background solution. The leading 
order derivative expansion for the effective action is obtained by first taking the 
(exact) expression for the effective Lagrangian in a covariantly constant background, 
substituting the space-time dependent background, and then integrating over space- 
time. For an instanton background, which is self-dual, we base our derivative expansion 
approximation on a covariantly constant and self-dual background |31I20||. This leads 
to the following simple integral representation for the leading derivative expansion 
approximation to the effective action ^Hj 

TiJA;m) = / — n41n 1 + — - + 7V3 — arctan^-— 

-84 + 7684 f 1, 1, ^; | - ^ Inm (24) 

Figure El shows a comparison of this leading derivative expansion expression with the 
exact numerical data. In the range covered, the agreement is surprisingly good for such 
a crude approximation. One would expect good agreement in the large m limit, but not 
in the small m limit. I conjecture that this remarkable agreement for both small and 
large mass is due to the special conformal symmetry of the single instanton background 
[31 El- This is physically sensible, but has not yet been proved. 

9. Concluding Comments 

Our angular momentum cutoff regularization and renormalization procedure of using 
the Gelfand-Yaglom result numerically for the low partial wave modes, and radial 
WKB for the exact large L behavior for the high partial wave modes, produces a finite 
renormahzed effective action which interpolates precisely between the previously known 
small and large mass limits. This is a highly nontrivial check, as all three computations 
(small mass large mass jHllHI, and all mass [Tj) are independent. Our result also 
provides a simple interpolation formula P] which can be translated into the instanton 
scale dependence at fixed quark mass, and this may be of phenomenological use. 

Our computational method is versatile and can be adapted to a large class of 
previously insoluble computations of one-loop functional determinants in nontrivial 
backgrounds in various dimensions of space-time, as long as the spectral problem of 
the given system can be reduced to that of partial waves. Hyunsoo Min's talk at this 
conference reports a new computation [23] of the false vacuum decay rate in a four 
dimensional self-interacting scalar field theory. Here the classical "bounce" solution is 
radial in 4d Euclidean space, so our technique is ideally suited to the computation of the 
fluctuation determinant about this bounce. An interesting feature here is the presence 




Figure 5. Plot of r-5(m), comparing the leading derivative expansion approximation 
(solid line) from (|24|l with the precise numerical answers (dots). The dashed lines show 
the small and large mass limits from ©. 



of a negative mode and zero modes. A different regularization method, also based on 
Gelfand-Yaglom, lias been applied to the false vacuum decay problem by Baacke and 
Lavrelashvili [221 ■ See also the very interesting earlier work of Isidori et al [221 ■ 

Finally, we note that our method provides an extension of the Gelfand-Yaglom 
result for ODE's to the case of separable PDE's. As noted by Forman ^3] for the 
2d radial disc problem, the naive extension does not work because the product over 
angular momenta diverges, even though the result for each angular momentum is finite. 
Physically, this is because in higher dimensions renormalization is required. Our method 
incorporates renormalization and yields a finite physical answer for the determinant of 
the separable partial differential operator. 
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